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Abstract 

Ph. The "normal" analysis of pp and pp elastic scattering uses a 'spinless' Coulomb am- 

Q_i| plitude, i.e., a Rutherford amplitude (2y/na/t) multiplied by a Coulomb form factor 

G 2 (t), an ansatz that pretends that the nucleon does not have any magnetic scattering. 
In this note, we investigate the role of the anomalous magnetic moment of the nucleon, 
k ~ 1.79. Given the method of analysis currently used by most published experiments, 
we conclude that the current experimentally inferred values of p for pp should be syste- 
meatically lowered by ~ 0.005 — 0.0100 and, correspondingly, the p values for pp should 
be systematically raised by the same amount. We discuss the theoretical uncertainties 
and a method of experimentally minimizing them. 



*Work partially supported by Department of Energy grant DA-AC02-76-Er02289 Task B. 



I Introduction 



We will only calculate electromagnetic amplitudes accurate to order a, i.e., one-photon 
exchange diagrams. Further, we will consider only high energy scattering (E^ ^> m, where 
m is the nucleon mass) in the region of small \t\, where t is the squared 4-momentum transfer. 
We will measure m and in GeV and t in (GeV) 2 , and will use h — c—1. 

La 'Spinless' Coulomb Scattering 

If we consider 'spinless' proton-antiproton Coulomb scattering, the relevant Feynman dia- 
gram is shown in Fig. ^| The electromagnetic differential cross section is readily shown to 
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Figure 1: One-Photon Feynman Diagram for 'spinless' Proton — Anti-Proton Coulomb Scattering, Pi+p^ — > 
Pf+p'f , with Couplings eG and — eG , where G{q 2 ) is the Electromagnetic Charge Form Factor of the Nucleon. 
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where the upper (lower) sign is for like (unlike) charges, t is the (negative) 4-momentum 
transfer squared, and m is the nucleon mass. 



For small angle scattering, the term (l — 
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At high energies, the correction term 
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over into the well-known Rutherford scattering formula, 
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1, so eq. (fj) goes 
(3) 



where the electromagnetic charge form factor G(t) is commonly parameterized by the dipole 
form 

1 



where A 2 = 0.71, if t is measured in (GeV) 2 . We note that this is the Coulomb amplitude that 
is normally used in the analysis of pp and pp elastic scattering, i.e., the 'spinless' analysis [p]]. 

I.b pp Scattering, Including Magnetic Scattering 

The relevant Feynman diagram is shown in Fig. |j, where magnetic scattering is explicitly 
taken into account via the anomalous magnetic moment k (~ 1.79). The fundamental 
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Figure 2: One-Photon Feynman Diagram for Proton — Anti- Proton Coulomb Scattering, Pi+p^ — » Pf+p'j, 
with Couplings eV^ and — eV^ , where V^ L — Fi-y^ + i^F 2 (r' J ^q vl with form factors Fi(q 2 ) and F 2 (q 2 ). 



electromagnetic interaction is 



K 



eV = e( + i—F 2 o^q u ) , </ = Pl - P , 



(5) 



which has two form factors Fi(q 2 ) and F 2 (q 2 ) that are normalized to 1 at q 2 = 0. The 
anomalous magnetic moment of the nucleons is k, and m is the nucleon mass. Because of 
the rapid form factor dependence on t, the annihilation diagram for pp scattering (or the 
exchange diagram for pp scattering) is negligible in the small \t\ region of interest and has been 
ignored. The interaction of eq. fl5|) is most simply treated by using Gordon decomposition 
and can be rewritten as 



eV = e 
Thus, using eq. (0), 
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the matrix element for the scattering is 
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where the upper (lower) sign is for pp (pp) scattering. A straightforward, albeit laborious 
calculation, gives a differential scattering cross section 
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We now introduce the electric and magnetic form factors, G_e(£) and Guif), defined as 

Gij(t) = F 1 (t) + —F 2 (t) and G M (*) = Fi(t) + kF 2 (*), 
and rewrite the differential cross section of eq. @ as 



(9) 



G?<7 

— = 47T 

oft 



a 



A 2 ab^ 

G|(t) 



x 



4m 2 



G\i{t) 



t 

iin 2 



1 + 



2m£- 



lab 



+ 



G\l(t) + 



1 (G 2 E (t)-GlM 



1 - 



4m 2 



+ G 2 M (t) 



8m 2 El h 



G%(t) 



t 

Am? 



G 2 M (t) t 



1 - sK 

4m 2 



2^ b 



(10) 



Note in eq. fllP]) that there is no cross term in GeGm, but only squares of these form factors. 
We can parameterize these new form factors with 
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where A =0.71, 
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with t in GeV/c 2 , and where G{t) is the dipole form factor already defined in eq. (f|), i.e., 
the form factor that is traditionally used in experimental analyses Q. 
We now expand eq. ( |10"D for very small \t\, and find that 
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where the new term in t, compared to eq. (fj), is 



2^ 



1 + 3.861*1 



3.39 



and is due to the anomalous magnetic moment of the proton (antiproton). To get an estimate 
of its effect, we note that G A {t) 1 — 11.26 1£|, in our units where t is in GeV/c 2 . We note 
that the new term is not negligible in comparison to the squared form factor, reducing the 
form factor effect by about 35% if the energy E^ is large compared to m. In this limit, we 
find that independent of the energy E 2 ah , 
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and is to be compared with the 'spinless' Rutherford formula of eq. 
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II Effects on Experimental Analysis of Elastic Scatter- 
ing 

UA4/2 has recently made a precision measurement || of p-p scattering at ^/s = 541 GeV, at 
the SppS at CERN, in order to extract the p value for elastic scattering. We now reanalyze 
this experiment, taking into account the magnetic scattering. They constrain the total cross 
section by an independent measurement @] of (l+p 2 )<7t t = 63.3±1.5 mb. For their published 
p-value of 0.135±0.015, this implies that they fix the total cross section at cr to t — 62.17±1.5 
mb. The main purpose of the UA4/2 experiment was the measurement of the p value, defined 
^ P = %mf ft=o) where f n (t = 0) is the forward nuclear scattering amplitude. 

II. a Spinless Analysis Neglecting Magnetic Scattering 

The experimenters parameterized the nuclear slope amplitude as /n(|^|) = e~ b],t ^ 2 , 

and measured the nuclear slope parameter as b = 15.5 ± 0.2 GeV -2 . They fit pp elastic 
scattering data at ^fs = 541 GeV over the t-interval 0.00075 < \t\ < 0.12 GeV 2 . They used 
for the Coulomb amplitude the 'spinless' Rutherford amplitude, modified by a Coulomb 
phase factor ia<f)(i), i.e., f c (t) = 2 ^ L G 2 (t)e ia ^ t \ where the phase||] </>(£) is given by 





where 7 = 0.577. . . is Euler's constant, b is the slope parameter, A 2 [= 0.71 GeV 2 ] appears 
in the dipole fit to the proton's electromagnetic form factor, G(t). The upper sign is for pp 
and the lower sign for pp. Using these parameterizations, the differential elastic scattering 
cross section is 

^° f 1 1 2 47ra 2 4 a a tot „M(+\\r>i(+\„-W\/i 1 °tgt(l + P 2 ) -b\t\ 
^ = |/ c + /n| ^—^-G(t) + —r--(p + a(p{t))G(t)e 1 17 + — e 1 (15) 

We now introduce the parameter t , defined as the absolute value of t where the nuclear and 
Coulomb amplitudes have the same magnitude, i.e., t = ^ = 14 Q ^ - , when cr to t is in mb, 
and to is in GeV 2 . For cr to t = 62.17 mb, we find that t = 0.00115 GeV 2 . We can now rewrite 
the differential cross section as 

W\ = ife { G4(t) ^ + 2 ¥\ ip + a ^ G2 ® e ~ bm + + p 2 ^ e ~ m } ' ( 16 ) 

for t min < t < t maX ) which was the form used by the UA4/2 group to analyze their experi- 
mental data. They extracted the value p = 0.135 ± 0.015, with a statistical error of 0.007. 
We emphasize that their analysis, using eq. ([16]), neglected the effects of the anomalous 
magnetic moments of the nucleons. 



II. b 'Spinless' Analysis, Taking into Account the Magnetic Scat- 
tering 

For our small \t\ analysis, we approximate G A (t) as G 4 (t) ~ 1 — 2a\t\, where a = -A = 
5.6338. We further write G 2 (t)e~ b ^^ 2 ~ 1 — (a + |)|t|. However, if we take into account 
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the anomalous magnetic moments, we see from eq. (13) that we could have written f c (t) ~ 
2-y/7r jj|G 2 (t) {1 + 1.93|£|} e 10 ^, an energy-independent result. Literally, we have used the 
'spinless ' ansatz that the Coulomb amplitude is given by the square root of the Coulomb 
cross section. Expanding the above equation in \t\, incorporating the new factor of 1 + 1.93 \t\, 
we can rewrite in a concise form the 'correct' Coulomb amplitude as 

fM*2^(&{t)<M £ >, (17) 

where 

G 2 ff (t) « 1 - a'\t\ = 1 - (a - 1.93)|t| = 1 - 3.674|t|. (18) 

Thus, a' = 3.674, and we can now write a modified form for the cross section which mimics 
eq. fll6l) as 

for t min < |t| < t max , by using G e g(t) in place of G(t), i.e., by replacing a by a' and p by p'. 

In order to extract the new value of p' from eq. (|19|) without having to refit directly the 
experimental data, we will require that the integral of eq. (|l!|) be equal to the integral of 
eq. ([TBI), from t min to t max . This insures that we fit the measured events with both formulae, 
allowing us to solve for p', given the UA4/2 published value of p = 0.135. Thus, we require 
that 

ax da' /-tmax da 

, G?|t| it mi T, 

To obtain an approximate analytical solution, we can set equal the contributions of the terms 
(1 + p' 2 )e~ b ^ and (1 + p 2 )e _fc '*', and expand in \t\ to first order. We find the approximate 
solution for p', the corrected value of p, valid when we set the contribution of a<p(t) = 0, 

pl = (a / -a)to + p(l-(a + V2) I ^^) 
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p+(a -a) t + p— if t < p — — - (21) 

\ m^max/ t-min J / m^max/ ''minj 

We will later see numerically that neglecting the a<p(t) contribution is a reasonable approx- 
imation, since the average value of a<p{t) = —0.00036 ~ over the t-range in question. We 
find, using the UA4/2 extracted values described above, that Ap = p' — p = p' — 0.135 = 
—0.0114. An exact numerical solution of eq. (|20|) , where we expand to first order in \t\ 
yields the final answer ApuA4 — —0.0105. This shift in p is comparable to the quoted total 
error of ±0.015 and is larger than the statistical error of ±0.007. The solution in completely 
insensitive to the slope parameter b, as well as to a<p{t). It mainly depends on the value of 
p that's measured, and critically on the t interval that was measured. We note that these 
conclusions are in contrast with the analysis of this problem made by N. Buttimore||, which 
did not take into account the experimental t interval. 
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Since all high energy p values come from experiments that measure about the same range 
of t and similar b values, it is likely that all the p values for pp scattering need to be lowered 
by ~ 0.005 — 0.010 and that all pp values need to be raised by about he same amount. For 
example, E710 measured§ p Pp = 0.134 ± 0.069, cr tot = 72.2 ± 2.7 mb and b = 16.72 GeV~ 2 , 
with t min = 0.00075 GeV 2 and t max = 0.077 GeV 2 , at ^ = 1.8 TeV. In this case, the exact 
numerical solution of eq. (20) gives ApE7io — —0.004. Thus, any global fit to the p values of 
pp and pp scattering should take these systematic shifts in the p values (at all energies) into 
account. 



II. c A Nuclear Model Taking Account of Magnetic Scattering 



In the preceding Section, we emulated the 'spinless' analysis by ignoring the structure of the 
nuclear scattering, and not distinguishing between spin flip and non-spin flip of the nuclear 
scattering, which we will now remedy. In this Section, we will introduce a toy nuclear model, 
which mimics the known electromagnetic matrix element of eq. (0), and adds coherently with 
it. We thus write the total matrix element for coherent nuclear and Coulomb scattering, 
taking into account the Bethe phase as expressed by CahnH, as 
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where now the upper sign is for pp and the lower sign is for pp scattering, since we introduced 
\t\ into eq. (B3). We have substituted e 2 = Ana and defined the strong coupling analog of 
a as h 2 = Ang, where h is the (complex) nuclear 'charge', kn is the nuclear 'anomalous 
magnetic moment' and H\(t) and H 2 (t) are the nuclear form factors. We have replaced the 
electromagnetic propagator |t| by the nuclear propagator m 2 . Later, we will use the optical 
theorem to fix the real and imaginary portions of h 2 and will also introduce "electric" and 
"magnetic" nuclear form factors Ge n and Gm n , analogous to eq. (0), with 
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The squaring of M[ in eq. (|22| ) will give rise to three terms, 
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since e %a ^ m 1 + ia<f)(t). The first term of eq. ( p5|) corresponds to pure Coulomb scattering, 
the last term to pure nuclear scattering and the term 2\M c \{^teM^ + a<p) to the coherent 
interference cross section between nuclear and Coulomb amplitudes. We note that the Dirac 
structure of all three terms in eq. ( f25|) is the same, which greatly simplifies the evaluation. We 
have already calculated the term \M C \ 2 and the substitution of a — > g, F\ — > Hx,F 2 — > H 2 and 
■Hi — > ^2 in the Coulomb term gives us the nuclear term |Mjv| 2 . Thus we find by inspection 



of eq. (|T2|) that the nuclear differential cross section is given by 
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We now expand eq. (|26|) for very small \t\, using eq. fl24]) and find that 
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Using the optical theorem, we now rewrite 
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since Gtv(O) = 1. Inspection of eq. 
value of the nuclear coupling is — — 
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differential scattering cross section of eq 



yields the equivalent statement that the (complex) 
(p + i) 9 ^. Using eq. ([28[), we can rewrite the nuclear 

as 
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The linear term in the brackets of eq. fl2"T|) or eq. ( ^9]) is clearly the spin flip term induced 
by the nuclear "anomalous magnetic moment" and goes to zero in the forward direction, 
as must be true for spin flip amplitudes. At ultra-high energies, eq. ( p?|) goes over to 
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Again, with a lengthy, but straightforward calculation, we find that the interference cross 
section is given by 
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Taking the limit of eq. (^) for small t, we find that 
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Introducing again the parameter t = we write, in the small t limit, the elastic differential 
scattering cross section for coherent Coulomb and nuclear interactions as 



da 



a 



a i 



1 + 



'k(k + 2) 



lab 

T 2^G\t)e-^ 

\t\ Plab 



2m 2 



2mE ] 



+ 



lab 



2£iab 



1*1 



p + a<j>(t) 



K + K N + KKat 



+ (l + p')e 



2 m 2 2mE\ a b 
( k n (k n + 2) 1 



1 + 1)' 



2-Efab 



1*1 



\ 2m 2 2m£i a b 
Finally, in the high energy limit, eq. (|3^) simplifies to 
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In order to use eq. (|3~3"D, we must expand the nuclear term for small t 



' da 



N 



Otot 
167T 

167T 

O-tot 2 
167T 



(1+P 2 



«jy(/Cjy + 2) U| 

+ 2m 2 111 
( _ k n (k n + 2) \ 
I 2m 2 J 



1*1 



(1 + p 2 )e~ b ' w , where 6' = b 



kn(kn + 2) 
2m 2 



(34) 



- 8 - 



Effectively, it is b' that is measured, not b. Thus we rewrite eq. (^) as 
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2m 2 



4m 2 
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We really have no deep knowledge of a value that is appropriate for k^, let alone its sign. 
Since the known polarization at Fermilab energies @, § is small, of the order of several per- 
cent, a magnitude compatible with it being produced by Coulomb interactions, this suggests 
that |/cjv| < k. It is tempting, however, to set k>n = k, since the conventional interpretation 
of the anomalous magnetic moment k is that it arises from the strong interactions. If we 
equate k and k^, we find the same interference term of eq. (|i~9"D where we "ignored" spin, 
i.e., 2G 2 (t)^e- b l'l/ 2 (l + 1.93|t|) . Thus, the numerical evaluat ion made earlier in Section lll.al 
is valid — namely, for UA4/2, the p value changes by « —0.01. 

On the other hand, even if we set kn = 0, we still would have an interference term 
proportional to (1 + 1.02|t|). 

We see from eq. (21) that for kn = that the shift in p, which is proportional to a' — a, is 
now reduced by 1.02/1.93 = 0.53. Thus, the p value for = is shifted by ~ —0.006. We 
conclude that for reasonable values of the parameter kjv, there is a significant shift in the 
p value of UA4/2 due to the anomalous magnetic moment of the nucleon. The theoretical 
uncertainty in the p-value can be reduced if the experimenters break their data up into two 
distinct regions — Region 1 being the interference region from t min to ~ 10 x t and Region 
2 from 10 x t to t max . In the case of UA4/2, this would change the uncertainty in ApuA4 
from « 0.006 to ~ 0.001, if we found p from Region 1 and cr to t and b from Region 2, even 
assuming that kn = 0. 

The other theoretical uncertainty is in the value of a$(£). The original derivation^ 
assumed that we had 'spinless' scattering, that the nuclear amplitude was e bt ^ 2 , where b' 
was the measured slope, and that the Coulomb amplitude was jG 2 (t). We see that we really 
should be using f G 2 ff (i), or, using A' 2 w 1.0 rather than A 2 = 0.71, in eq. fli~i]) . Fortunately, 
this is a very small change, and increases the pp p-value by ~ 0.001, thus contributing 
negligibly to Ap UA4 . 

In conclusion, it seems sensible for experimenters to redo their data analysis using = k, 
i.e., using G e g(t) and A' 2 rather than G(t) and A 2 , in the two distinct t-regions described 
above. This procedure allows the experimenter to control the theoretical uncertainties, i.e., 
our lack of knowledge of the nuclear amplitudes. 
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